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Elementary excitations of incompressible quantum liquids (IQL's) are anyons, i.e., quasiparticles 
carrying fractional charges and obeying fractional statistics. To find out how the properties of these 
exotic quasiparticles manifest themselves in the optical spectra, we have developed the anyon exciton 
model (AEM) and compared the results with the finite-size data for excitons of i/ = 1/3 and p — 2/3 
IQL's. The model considers an exciton as a neutral composite consisting of three quasielectrons and 
a single hole. The AEM works well when the separation between electron and hole confinement 
planes, h, obeys the condition h > 21, where I is the magnetic length. In the framework of the AEM 
an exciton possesses momentum k and two internal quantum numbers, one of which can be chosen 
as the angular momentum, L, of the k = state. Charge fractionalization manifests itself in striking 
differences between the properties of anyon excitons and ordinary magnetoexcitons. Existence of 
the internal degrees of freedom results in the multiple branch energy spectrum, crater-like electron 
density shape and 120° density correlations for k = excitons, and the splitting of the electron 
shell into bunches for k 7^ excitons. For h ^ 21 the bottom states obey the superselection rule 
L = 3m, m > 2 are integers, and all of them are hard core states. For h ~ 21 there is one-to-one 
correspondence between the low-energy spectra found for the AEM and the many-electron exciton 
spectra of the v = 2/3 IQL, whereas some states are absent from the many-electron spectra of the 
u — 1/3 IQL. We argue that this striking difference in the spectra originates from the different 
populational statistics of the quasielectrons of charge conjugate IQL's and show that the proper 
account of the statistical requirements eliminates excessive states from the spectrum. Apparently, 
this phenomenon is the first manifestation of the exclusion statistics in the anyon bound states. 
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I. INTRODUCTION 

Incompressible quantum liquids (IQL's)0 underlie the fractional quantum Hall effect (FQHE) discovered by Tsui, 
Stormer, and Gossarda The chaijee carriers in these liquids are anyaps, i.e., quasiparticles (quasielectrons and quasi- 
holes) carrying fractional chargea^l and obeying fractional statistics.QQ Historically the main experimental discoveries 
in this field were done by magnetotransport experiments. However, the role of spectroscopic methods is continuously 
increasing since they provide an indispensable tool for investigating spectra of elementary excitations. Fine structure 
specific for different electronic phases was discovered in the spectra of radiative photoemission.B These spectra were 
used to measure gaps in the energy spectra of IQL's, for investigation of phase transitions between the IQL and 
Wigner solid phases, etc. The frequency of long-wave neutral elementary excitations of IQL's was measured in Raman 
scattering experiments.l3 

A challenging problem in physics of IQL's is a direct observation of the charge fractionalization. Between different 
exciting approaches to this problem the spectroscopic approach seems to be one of the most promising. Indeed, 
spectroscopy permits one to observe properties of IQL's in the bulk where the effect of the impurities and edges 
is reduced to the minimum. The intrinsic spectroscopy of IQL's is the magnetospectroscopy of excitons. However, 
properties of excitons reflect the spectrum of the elementary excitations of an IQL (quasielectrons and quasiholes, 
magnetorotonsjj etc.) and can be treated in terms of them only when the separation h between electron and hole 
confinement planes is sufficiently large. Indeed, when h I, where / is the magnetic length, the filling factor v of 
the liquid strongly deviates in the vicinity of the hole from the quantized value v = p/q. Under these conditions 
the properties of the IQL cannot be treated in terms of its quasiparticles. The spectroscopy of a remote hole has 
been discussed from different standpoints in a number of papers Jail3 Despite the fact that experiments in the h ^ I 
region are rather complicated, the separations up to ft. « bl were achieved in experiments on extrinsic radiative 
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photoemission.Ei Of special importance might be ^experiments performed for a fixed filling factor v = p/ q and a 
variable dimensionless electron-hole separation /i/Z£j The first experiments of this kind were reported recently.LJ 

The investigation of excitons is also important from the different standpoint. In the theory of IQL's the statistical 
properties of the system of free anyons are usually discussed. In an exciton the anyons exist in a bound state because 
of the attracting potential of a hole. It was showrJlSllII that the energy spectra of excitons of = 1/3 and v = 2/^ 
IQL's are closely related to the statistical and dynamical properties of quasielectrons of these IQL's. Therefore, the 
exciton problem can be really treated as a few-anyon problem. Excitons represent a wide class of few-anyon systems. 
It includes excitons, anyon ionsEErO (which can be also treated as quasiparticle-exciton complexescj) , some impurity 
centers, etc. In what follows, we consider excitons as bound few-anyon composites and investigate their properties in 
some detail. 

There are two approaches to the theory of excitons in IQL's. The first approach is based on the exact diagonalization 
for finite-size systems in the spherical geometryEI Having in mind workable system sizes, it provides reliable results 
for h ^ 2Z, at least as applied to the v = 1/3 and-ii-j^ 2/3 IQL's. The second approach is based on the anyon exciton 
model (AEM) proposed by the present authorsJlil'E3 In the framework of this approach an exciton is considered as a 
neutral composite quasiparticle consisting of several anyons and a hole. This model is exact only when the exciton 
size which is about h is large compared to the size of anyons which is about /, i.e., for h/l 3> 1. Therefore, the two 
approaches are complementary and one can expect that they match when h/l 1. 

It is the main statement of the AEM that excitons of IQL's possess a multiple branch spectrum. Indeed, a charged 
particle at the lowest Landau level possesses a single degree of freedom and a single quantum number. Aa exciton being 
a neutral entity possesses in a magnetic field a vector momentum k absorbing two degrees of freedom.cJ Therefore, an 
exciton consisting of q anyons and a hole possesses g — 1 internal degrees of freedom. For an ordinary magnetoexciton 
q = 1 and the spectrum consists of a single branch. For q > 2 an exciton acquirfts-|iBtp»ial quantum number(s) 
and multiple branch spectrum. This prediction of the AEM permitted Apalkov et aZ.llaO'cil to represent the energy 
spectra found by finite-size computations for = 1/3 and = 2/3 IQL's as a system of exciton branches and to assign 
to these branches internal quantum numbers. The latter detepmine the values of the exciton angular momentum, L, 
in the k = states. Zang and BirmarEj and Chen and QuinncB also inferred existence of several exciton branches in 
their finite-size data. 

In the range accessible for finite-size studies, h ^ 21, there are two types of excitons in the lower part of the energy 
spectrum, anyon excitons and tight excitons. Anyon excitons are loose entities with a pronounced anyon-deasity dip 
at the center. They are generically related to the quantum states making up the low-energy (anyon) sectortZlEEl of the 
electron subsystem. Each anyon exciton is a bound state of a three-quasielectron complex from the anyon sector and 
a hole. This finding establishes a connection between the spectroscopy of excitons and the low-energy physics of the 
FQHE. There is a striking difference in the quantum numbers of the low-energy anyon excitons of the = 1/3 and 
J' = 2/3 IQL's. It is related to the difference in the energy spectra of the three-quasielectron complexes originating 
due to the different populational statistics of the quasielectrons of these liquids. Tight excitons are dense entities. A 
sharp density maximum is achieved either in the center of an exciton or in a close vicinity of the center. These excitons 
are not related to the low-energy sector of the electron subsystem. The L = tight exciton is the bottom state of the 
exciton spectrum in all the region h 21. Therefore, the h ^ 21 region seems to be of the most importance for the 
study of the anyon substructure of excitons and for the spectroscopic observation of the charge fractionalization. 

We suppose everywhere in this paper that the spin-polarized background is stable with respect to formation of 
charged spin-textures (skyrmions). It was predicted recentlyo that spin-textures develop in the ground state of a 
1/ = 1 exciton for large values of h/l; the critical value of h/l is about 1.5 for GaAs. Nevertheless, we feel that 
the spin-polarized ground state-iof a, v = 1/3 magnetoexciton is stable because skyrmions can exist only at very low 
magnetic fields when v w 1/3.E3 Spin-depolarized excitons should be seen in the spectrum of excited states. 

In this paper we develop the theory of excitons in the framework of the AEM as applied to excitons consisting 
of three quasielectrons and a hole. We use Halperin pseudo-wave-functionsH in a boson basis and investigate the 
energy spectruffi|,[-e|lectron density distribution, and anyon correlations in an exciton. We make comparison with the 
finite-size datatjO on the excitons of the v = 1/3 and v ~ 2/3 IQL's and conclude that excitons of the v — 2/3 
IQL are described rather well by the AEM because of the bosonic populational statistics and the narrow form-factors 
of the quasielectrons of this liquid. We also relate the difference in the energy spectra of the anyon excitons of the 
I' = 1/3 and = 2/3 IQL's to some specific features in the exciton shape found in the framework of the AEM. For 
both IQL's the bottom exciton states are made by tight L = excitons for h < 21 and by a succession of anyon 
excitons for h ^ 21. This succession consists of hard core excitons with the angular momenta L which are integers of 
3 and increase as L oc /i^ with h. 

The general outline of the paper is as follows. In Sec. |^ we construct a full basis of the exciton wave functions for 
the AEM using the translational symmetry and permutation symmetry arguments. In Secjn wc develop a technique 
for calculating different matrix elements entering the Schroedinger equation. In Sees. IV, M, and VI we obtain energy 



spectra, electron density distributions, and the density correlation functions, respectively. The latter functions unveil 
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the anyon substructure of excitons both in the finite-size data and in the AEM. In Se c.|VIl| we make a comparison of 
the results obtained in the framework of the AEM with finite-size data of Refs. ^s] and |17|. We propose that a striking 
difference in the finite-size data for the v — Xji and = 2/3 IQL's originates from the difference in the populational 
statistics of the quasielectrons of these liquids. 



II. WAVE FUNCTIONS 



Let us consider an exciton consisting of a valence hole with a charge (-l-e) and three QE's with electrical charges 
(— e/3) and statistical charges a. Such an entity provides the AEM description of the anyon excitons of the v = 1/3 
and V = 2/3 IQL's. For q.vl v — 1/3 IQL the statistical charge equals a — ^ l/i^i while for a,n v = 2/3 IQL the 
statistical charge has the same value, a = 1/3, as for quasiholes in an i/ = 1/3 IQL.U In comparison, a = for bosons 
and a = 1 for fermions. In the strong magnetic field limit, when the Coulomb energy ec = e'^/el <^ fiuJc, where ojc is 
the cyclotron frequency and e is the dielectric constant, it is convenient to employ dimensionless variables scaled in 
units ec, I and e. We use the symmetric gauge, A = z x r/2, where z is a unit vector perpendicular to the confinement 
plane. Instead of the hole, r^j, and anyon, r^, coordinates it is convenient to introduce the following two-dimensional 
(2D) coordinates: 

1 1 . ^ 1 

= 2 + 3 X! ''^^ ' ^^^j ^ ^ji = ~ = 23, and 31 . (1) 

i=l j=l 

R has a meaning of the center-of-mass of the exciton coordinate. The coordinates p and rji are the internal variables 
which are not affected by the translational motion of an exciton. Complex coordinates Zji = Xji -\- iyji, as well as rjj, 
are not independent. Indeed, 

ri2 + r23 + rsi = . (2) 

Despite the fact that the constraint ^ results in some complications, the introduction of the variables Yji enables 
one to develop the theory in a form symmetric in all anyons and, therefore, finally simplifies the equations. 

Anyons and hole live in two different parallel planes separated by the distance h. Nevertheless, only the 2D 
coordinates of Eq. |l| enter into the excito n wav e function. The separation h enters only into the Hamiltonian of the 



anyon-hole interaction derived in the Sec. Ill C 



The most general form of the pseudo-wave-function of an anyon exciton meeting all general requirements is as 
follows: 

^-L.kCR, P, {%}) = cxp{ikR + Iz • (p X R) - l(p - df}PL{. ■ ■ zji . . .) 

X n (^ji)"exp{-|z,j|V36}/%/2^ , (3) 

where the pair of indices, jl, takes the values specified in Eq. (|l|), A is the normalization area, and Pl is a homogeneous 
polynomial in coordinates Zji of the degree L. 

The basic properties of the functioiis-|5'L,k can be checked by inspection, (i) Since the exciton is a neutral entity, 
it possesses an in-plane momentum k,c3 and ^'i.k satisfies the equation of magnetic translations: 

Ta^-L.klR, P, {zji}) = exp[i a • A(p)]«'i,k(R - a, p, {z,,}) = e-''' =^^'i,k(R, P, {%}) • (4) 

The parameter d = z x k is related to the dipole moment of the exciton, (— d). (ii) The function VP^.k belongs to the 
lowest Landau level. Indeed, the nonanalytic factor of it can be shown to have the form exp{— J^j ~ il-^/tP}, 
whereas the other factors are analytic functions of Zj and Zh- (iii) The function Pl{. ■ ■ zji . . .) is a homogeneous 
polynomial of the degree L which is symmetric in all coordinates Zj. These polvjiomials form a boson basis, and 
the effect of the fractional statistics is taken into account by the factor J|^;(zj;)".aEil (iv) A system of four charged 
particles in a magnetic field possesses four quantum numbers. Two of them are absorbed in the 2D momentum k. 
Two others determine the form of the polynomial P^ and are internal quantum numbers of an anyon exciton. The 
operator of the z projection of the angular momentum commutes with the Hamiltonian and the square of the 
momentum, k^, but it does not commute with the projections of k, i.e., with kx and ky. Therefore, the function '^l.Is. 
chosen in the k^, ky representation is simultaneously an eigenfunction of only for k = 0, and in this limit Lz — —L. 
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Therefore, the quantum numbers of an anyon exciton include the 2D momentum k and the projection of the 
angular momentum, = —L, of the exciton with k = 0. The angular momentum L numerates branches of the 
exciton spectrum. The fourth quantum number, which will be specified in what follows, numerates branches with 
coinciding values of L. The multiple branch structure of the anyon exciton spectrum is a direct consequence of the 
charge fractionalization which results in the appearance of the internal degrees of freedom of an anyon exciton and of 
the related internal quantum numbers. 

By definition, the polynomial Pl is symmetric in coordinates Zj . To establish the symmetry of it in the symmetric 
coordinates z^/, one can start with a monomial z'2^23^311 ^-Pply to it all operations of the permutation group, and 
take the sum over the group. This transformation results in the polynomial 

( ~h 5^2 I 5^3 5^2 1 5^3 5^1 \ 1 ( 1 +'2+^3 (' 5:'3 5^1 1 5^1 5^2 i 5^3 ~h :ph \ 

1^12'^23-^31 "1" -^23-^31^12 -'31^12-'23-' I ) ^, -'12-^31-^23 -'23^12^31 ^31-^23^127' 

which has different properties depending on the parity of L = + Z2 + ^3- If i is even, the polynomial is a permanent, 
and, therefore, is symmetric in the coordinates Zji. However, when L is odd, the polynomial is a determinant which 
is obviously antisymmetric in the coordinates Zji and is non-equal to zero only for h h 7^ h- For example, for the 
lowest possible value oi L, L = 3, this determinant turns into a Vandermonde determinant 



W{zi2,Z23,Z3i) 
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?2 


^12 


-*23 


^31 



(Zi2 - 2:23)(-Z23 " ^3l)(2;31 " ^12) • (5) 



Therefore, L-even and L-odd polynomials have rather different properties. All of them are symmetric in bosonic 
permutations zi <-> Z2, etc., but they have opposite symmetry with respect to the permutations of the zi2 Z23 type. 
To find the explicit form of the polynomials Pl , it is convenient to introduce new real coordinates 

1 ^ 

^3 = rj - ro , ro = - ^ r^; , (6) 

1=1 

where rp is the center-of-mass of the anyon sybsystem. The corresponding complex coordinates are Q — £,jx + i^.jy 
These coordinates are subject of the constraint: 

3 3 

E^. = 0' Eo = o. (7) 

In these coordinates the bosonic symrnetry of Pl has the usual form, and we can apply the fundamental theorem 
of the theory of symmetric polynomials.E3 According to it, Pl {(,11^2X3) can be expressed in the unique way as a 
polynomial in the three elementary symmetric polynomials: 



Pl = Cl + C2 + C3 , P2 = C1C2 + C2C3 + C3C1 1 P3 = C1C2C3 ■ (8) 

The first polynomial is equal to zero, pi = 0, because of the constraint of Eq. (^). Therefore, the polynomials Pl are 
polynomials only in p2 and ps. The first L-even polynomials are: Pq = const, P2 cx p2, Pa oc p\^ i.e., there exists only 
a single elementary polynomial of a given degree L. However, two elementary polynomials, p\ and p|, contribute to 
Pg- It is easy to check that the number of basis functions increases by one each time when L takes values L = 6m, 
where m is an integer. Therefore, the number of L-even polynomials is equal to [L/6] + 1, where [L/6] is the integral 
part of L/6. All L-odd polynomials can be obtained by multiplying L-even polynomials by p^. The latter equals 
P3 = —^W because of Eq. (^). 

We are now in position to choose a full basis of polynomials Pl in symmetric coordinates Zji. These coordinates 
are most convenient for all the following calculations. L-even polynomials can be chosen as 

p _ -L-4M -2M -2M , -L-4M -2M -2M , -L-4M -2M -2M fQ\ 
^L,AI — Z12 ^23 ^31 + ^23 ^31 ^12 + ^31 ^12 ^23 J Wl 

where M — 0, 1, . . . [L/6]. Polynomials Pl.m are linearly independent, and the total number of polynomials with a 
given L is equal to [L/6] + 1. All linearly independent L-odd polynomials can be obtained as: 

Pl,m = WPl-sm , ^3,0 = W . (10) 
The total number of them equals [(L — 3)/6] + 1. 
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This choice of polynomials determines the full set of quantum numbers in the wave function of Eq. (||) as L—M, and 
k. To our best knowledge, in the previous studies only the L-even polynomials have been taken into accountS When 
choosing polynomials Pl,m, we have not imposed the hard core constraint and defer the discussion of the related 
properties to what follows. 

It is an important feature of the AEM that the wave functions (and, therefore, electron densities, etc.) of all 
eigenstates with L < 5 and also L = 7 arc completely determined by the symmetry requirements. They do not 
depend on the specific form of the Hamiltonian and, in particular, on h. 

III. THE SCHROEDINGER EQUATION 

In this section we calculate the Hamiltonian of the AEM in the Vl/L,M,k basis as the matrix of a point charge Coulomb 
interaction. It means, that we aa^ect form-factors of quasielectrons which have a scale of several magnetic lengths and 
are known only approximately.EjE3 We postpone the discussion of inaccuracy originating from this approximation to 



Sec. VII 



Unfortunately, functions 'i^ L,M,k are orthogonal only in quantum numbers L and k. As a result, the scalar products 
(^L,M,k, ^'L,Af',k) 7^ for M 7^ M', and the matrix B of these scalar products is block-diagonal. The size of blocks 
is equal to 1 for L < 6 and L — 7 and increases by 1 each time when L increases by 6. With a non-diagonal matrix 
B the Schroedinger equation has a form: 

Hx^eBx, (11) 

and one has to find matrices H and B. To perform the calculations, it is convenient to employ variables R, p, and 
three r,jj and to take into account the constraint of Eq. (|^) by the usual transformation: 



(5(ri2 + r23 + rai) = j exp{if • (ri2 + r23 + 



)} • (12) 



It adds the new variable f , but all calculations become symmetric in anyon variables. The Jacobian of the transfor- 
mation is equal to 1. 

The Hamiltonian is diagonal in k, therefore, we write out only the diagonal in k matrix elements. For k = 0, the 
Hamiltonian is also diagonal in L, and Eq. (O) acquires a block-diagonal form. 



Since all terms in the polynomials Pl,m,\l , Eqs. (H) and (10), have the same form, we concentrate in what follows 



on the matrix elements taken in the basis of the functions 5'{„}k : 
*{n}k(R, P, {%}) 

= exp{zkR +^z-{pxR)- i(p - d)2}zt|+"z2"3i+"z3",^+" exp{- ^ \z,i\y36}/V2^ . (13) 

jl 

Here polynomials Pl,m are substituted by monomials, and {n} denotes a set of quantum numbers rti,rt2 and 713. In 
the next parts of this Section we describe in some detail the technique for performing different types of integrals. 

A. Nonorthogonality matrix elements 

The scalar products of functions ^'{„}k, when written in the variables R, p, and r^/, have the form: 

B{n}{n'} = ({"} I W}) 

^ I I '^'^ I (^^ / ^''I2(^'"23dr3i4'{„}k*{n'}k cxp{if • (ri2 + r23 + r32)} . (14) 

Integration over R and the Gaussian integration over p are straightforward. Since ^'{„}k is multiplicative in the 
variables r^/, one can rewrite Eq. ( p^ ) in the form: 

3 

(wiM)-/ (|^n<^(f). (15) 
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where 



M^l,{i) ^ j dr e-'-'/i8+«f-r (16) 

Since 

dip e±*'"'^+*^''™"'^ = 27ril™l J|„|(/r) , (17) 

where J\m\{fr) is a Bessel function, the angle integration in ( p^ results in: 

^il'(f) = 27rzl™-™'l exp[z(m - m')^f] A^Ll' W , (18) 
where i = 9/^/2, (^f is the azimuth of f, and 

M^±,{t) = / dr ri+2"+™+'"'e-'-Vi8j|„_„,,|(/,) . (19) 
Jo 

This integral can be expressed in terms of the confluent hypergeometric function <i>(/3,7;t) as0 

n\rn-\-m |/2+a o |m+7n' | +2(a+l) 

■^Ll'(i) =r(max{m,m'} + a + l) \ ^ ^ \i 

X $(niax{m, m'} + a + 1, \ m - mf \ + 1; -t) . (20) 

Here niax{m, m'} is the larger of the integers m and m' . After the integration over ipf in Eq. (|l5|), the coefficients 
({rt} I {n'}) take the form: 



({n} I {n'}) = J, 



A'^^/i? rfin*'"^""^'-^si'W ' (21) 



where 

n = ni + 712 + 713 • (22) 

Therefore, the scalar product ({71} | {n'}) of two functions ^{„}k is reduced to an one-fold integral from the product 
of three confluent hypergeometric functions. Matrix elements ({71} | {71'}) do not depend on k. Scalar products 
(L, M I L , M ) of two functions L,M,\i , which include polynomials Pl,m, are linear combinations of the coefficients 
{in} I {71'}). They do not depend on k, too. The coefficients of these combinations can be found from the equations 
(|), (|), and (0), but the final expressions are rather cumbersome, especially for L-odd polynomials. Therefore, we 
do not write out here their explicit form. r— ■ 
Integrals ( ^l| ) can be simplified for a ~ 0. Indeed, in this case the Kummer transformation&Zl 

$(/?,7;0 = e*$(7"/3,7;"0 (23) 

results in a <& function with the first parameter (7 — f3) equal to a negative integer. This function reduces to a 
polynomial, and $(/3,7;i) to a polynomial multiplied by e*. Therefore, the integral in Eq. ( pT]) for {{n} \ {n'}) can 
be performed exactly. This transformation highly simplifies calculations. Indeed, in the large h region, where the 
criterion of the applicability of the AEM is satisfied, the statistical parameter a can be neglected as it is shown in 
Sec. below. 



B. Anyon-anyon interaction 



The Hamiltonian of the anyon-anyon interaction is: 



Vaa^l{\zi2\-' + \z23\-' + \z3l\-'}- (24) 
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Matrix elements of Vaa in the basis of the functions ^{n}]i can be calculated by analogy with the matrix elements 
{{n} I {n'}). The denominators \zji\ ^ lower the power of r by 1 in one of the (f) factors entering in Eq. dlj). 

The final expression for the matrix element is: 

({n} I I {n'}) = i J ^{Mi:„7)(f)Af(i, (f)M(i, (f) 

+K:k (f (f (f ) + (f )<«7^ (f )} • (25) 

Like the nonorthogonality matrix elements, matrix elements of Vaa also do not depend on the momentum k. For this 
reason, the matrix of the operator Vaa is diagonal in the angular momentum, L, for arbitrary values of the momentum 
k. 



C. Anyon-hole interaction 

The Hamiltonian of the anyon-hole interaction has a form: 

1 ^ 

Vah = - :t E ^^-'^ ' ^^■''(''^■'^) = "^Jh > (26) 
where the three-dimensional anyone-hole separation rih should be expressed in terms of the difference coordinates: 

ri/i = p+^(ri2-r3i) + z/i . (27) 

Similar equations hold for r2h and r^h- 

It is convenient to introduce the Fourier image 

Vihirih) = j (^^^°''(9) exp{«q • P + • (1-12 - rai)} , (28) 

where Vahil) — (27r/g) exp(— g/i). The integrations over R and p in the matrix elements of Viui'Tih) can be performed 
in the same way as in Eq. (|l4|), and the result can be expressed in terms of the coefficients M^^, , Eq. (|l^), as: 



(M I v.. I {n'}) ^--Jj^Jj^ VM c 



X <U (f ) (f - <»3 (f + ^/3) ■ (29) 

After the angular integration, the factors M^'^\ take the form of Eq. (^8|), and Eq. (^) can be rewritten as: 

exp{j(/jf (ni — n[) + i(p^{n2 — n'2) + i(p+{n3 — rig)} 



■^l?„;(9/72) MZ'jm/^) M'Z'^{%fj2) . (30) 



I*-"^ (9/^/2) A^*-"^ ^Of2/o^ (-0/2 

Vectors f± are defined by the equation 



f±=f±q/3. (31) 
One can perform one angular integration in Eq. (^0|) if the phase (^f is eliminated by change in the variables: 

= <y9f - yjq , tp± = ^P±- Iff . (32) 

Here </?± are phases of vectors f± . 

After the integration over ^pc^^ Eq. ( p9| ) takes its final form: 
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1 nOO poo 

({«} I Vih I {n'}) = -- rr / dff / g e-''/V|„_„,|(fcg) 

•J j-^i Jo Jo 

X dip exp{i(^(n — n') + i'0_(n2 — ^2) + j'0+(?i3 — ^3)} 

■/o 

X ■^it;(9/V2) A^it. (9/^/2) A<il. (9/V2) ■ (33) 

Here the direction of k was chosen along the x axis to ehminate the complex phase from the matrix element (|3^). 
Explicit expressions of and f± in terms of the integration variables are as follows: 

e^'^± =(/±|e-^^)//± , /| = /2 + (g/3)2±^/gcos^. (34) 

One can check by inspection that Eq. ( [33[ ) is symmetric in the indices '^Si'^s: ^^^id that matrix elements 

are real for all sets of the quantum numbers {rtj}, compatible with our choice of the polynomials Pl.m, Eqs- (||) 

and ([l0|). Matrix elements of the operators V2/1 and V3/1 can be written by analogy with Eq. (^3|). 

Equations ( p5| ) and ( ^3|) for the Hamiltonian H and Eqs. ( |l4[ ) and ( ^l|) for the matrix i? determine completely the 
Schroedinger equation (pJj). 



IV. ENERGY SPECTRUM 



For k = the operator of the anyon-holc interaction Vah becomes diagonal in the angular momentum L. Therefore, 
the Schroedinger equation ( pj] ) is also diagonal in L since Vaa and B are diagonal in L for arbitrary k. It was shown 
in Sec. |l| that for each value of L such that L < 5 01 L — 7 there exists a single eigenfunction, and it does not depend 



on h. For these values of L the equations of Sec. Ill immediately give the energies of k = states. When the number 



of Pl,m polynomials with a given value of L becomes two or more, equations of Sec. HI give the coefficients of secular 
equations of the second, third, etc., order which determine the h dependent eigenfunctions and eigenvalues. 

The order in which k levels are arranged changes with increasing h. The main regularities can be understood 
using classical arguments. In the classical limit, which is achieved for a large exciton size, the exciton ground state 
takes the shape of an equilateral triangle with anyons in the vertices and a hole in the center. The anyon-anyon 
distance in this triangle, ri2, found from the minimum of the electrostatic energy, is equal to 

ri2 = ^/3/2h . (35) 

It will be shown below that triangular configurations are described by the polynomials Pqm,m and W Pqm,m- A 
straightforward calculation based on equations (|^) , (^ and ( p^ ) shows that a mean-square value of the inter-anyon 
distance in these states is equal to 

<rl^>=Q{L + 2 + ?,a) . (36) 

Comparing equations ( ^ ) and (|3^) one finds that 

ri2-/i-ii/2 (37) 

for L ^ 1. Therefore, with increasing h both the size, ri2, and the angular momentum, L, of the exciton ground 
state increase. This means that the order of the energy levels changes, and the bottom state possesses the angular 
momentum L ^ . Since for k = the Hamiltonian is diagonal in L, the level interchange occurs usually as a level 
crossing. 

In Fig. |l] the ground state energy is plotted as a function of h for two values of the statistical charge, a ~ —1/3 
and a = 0, by curves A and B, respectively. It is seen that both curves show the same gross features, including 
increase of the ground-state energy and the angular momentum with h. However, fine details are very different in 
the h <2 rcgjrm- We do not discuss these differences in more detail since they are expected to be sensitive to anyon 
form-factorsEJ^ES which were not taken into account in our calculations. However, since curves A and B practically 
coincide for ft, > 2, we believe that in this region the AEM provides reliable results. In what follows we restrict 
ourselves with this region and neglect the statistical charge, i.e., consider the bosonic model, a = 0. All data below 
are presented for this model. It was shown in Sec. Ill A, Eq. (|2^), that for a — matrix elements A^mm/(t) can be 
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expressed in terms of elementary functions. This fact permits one to reduce the three- fold integral of Eq. ( p3| ) to the 
one-fold integral: 

/"OO 

< LMk\Vah\L'M'k >= - / exp(-3gV2 - qh) J\L-L'\{kq) QLM,L'M'{q) dq , (38) 

Jo 

where functions QLM.L'M'{q), real and symmetric in indices, are polynomials in q. The lower polynomials are of a 
simple form: Qoo.oo = 1, Q20,oo = 9^/2, Q20,20 = I — q^ +(7''/4. Application of the equation ( |3^ ) highly simplifies all 
computations. 

For comparison, in Fig. |l| is also shown the energy Sd = — (2/3)'^/^/ft- calculated in the classical limit, Eq. (|35|). It 
is seen that in the region < 4 it differs considerably from the exact quantum data. 

Since for k 7^ the term Vah is non-diagonal in L, the dispersion law e(fc) can be found only numerically. In Fig. |^ 
it is shown for two values of h. The basis of polynomials used in computations included 22 L-even polynomials with 
L < 18, Eq. (^, and 22 L-odd polynomials with L < 21, Eq. (^. The following regularities are distinctly seen. As it 
was argued above, the levels with higher L values draw closer to the spectrum bottom with increasing h. The level 
interchange manifests itself as avoided level crossings. The level splitting near these crossings increases with k and 
decreases with the difference \L — L'\. These regularities can be understood if one takes into account that the Bessel 
function in the integrand of Eq. ( |3^ ) shows the power-law behavior, J\L-L'\{kq) oc {kq)^^~^ in the small k region. 
In particular, for \ L — L' \ = I the interaction of two branches can result in the negative exciton effective mass for 
small k values. Negative exciton dispersion can appear even in the ground state as it can be seen in |Eig. |^a. It is 
interesting to mention that exciton dispersion near k = is always positive in the two-semion problemJiil 

The above results show that the charge fractionalization determines both i) the basic multiple-branch structure of 
the exciton energy spectrum and ii) numerous specific features of the spectrum including the h dependence of the 
arrangement of the branches, avoided branch intersections, etc. 



V. ELECTRON DENSITY 



The distribution of the electron density, D\[r, k), around a hole can be found from the equation: 

1 ^ 

Z?A(r,k) = (-^J(r, -r,,-r)>Ak , (39) 

where the averaging is performed over the quantum state (A,k), where A numerates exciton branches. The density 
£'A(r,k) is exactly the quantity which i) permits one to check reliability of the model and in which ii) the specific 
pattern of the charge fractionalization manifest itself in a rather spectacular way. It is one of the basic criteria of the 
AEM that the excess charge density is small compared with the density of the IQL, v/2tt, for < 1/2 (and compared 
with (1 — v)/2tt for v > 1/2). Therefore, for v = 1/3 and 1/ = 2/3 IQL's one can expect that the AEM becomes 
applicable only when D\ <C I/Qtt ~ 0.05. This criterion will be applied in what follows. 

The explicit expression for D\ can be obtained in the same way as Eq. (^8|). Indeed, the operator Vah and the 
operator of the electron density of Eq. ( |39| ) depend on the same arguments, — Th- Therefore, the integrands differ 
only in the substitution Vah{q) by the Fourier image of S{rj — Vh). The final expression is: 

^^(^'1^) = ^ E ^os[{L-L')0] 

LM.L'M' 

X dqq exp(-3(zV2) J\L-L'\{q\r ~ d|) Qla/,L'A/' (<z) xi'Af (k)xiM(k) , (40) 
Jo 

where 6 is the angle between the vectors d — r and d. Dx{0, k) shows the electron density on the hole. 

For k = 0, the density distribution Dxir, 0) is shown in Fig. |^ for L < 6. In this case A can be completely identified 
by the index L for i < 5, but there are two functions for L — Q. For L < 5 the energies, the eigenfunctions, and 
therefore also the densities Dj^, do not depend on h. The L = state has a high density; Dq{0, 0) > 0.05. Therefore, 
the shape of the curve DqIt, 0) can not. bfi rfiliab le. Nevertheless, it is^emarkable that numerical calculations performed 
in the spherical geometry for v = 1 / 3L3llj'E3Ea and v = 2/3 IQL'alll convincingly show that for h ^ 2 the spectxiim 
bottom is made by the L = exciton having the electron density which is very close to the Fermi limit, l/27rJl3'EHl 
Functions Dl,m become broader with increasing L, and I?L,Af (0,0) decreases. The state L = 3 from which L-odd 
polynomials originate is especially remarkable because it is the first to show a crater-like density distribution with a 
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minimum at r = 0. This minimum is a signature of the charge fractionaHzation, since for ordinary magnetoexcitons 
the density shows a maximum at r = 0. We will discuss the properties of this state in more detail in Sec. VI. For 
L — 6 there are two eigenf unctions; they depend on h. In Fig. ^they are shown for h = 0; the lower energy component 
is drawn by a solid line. 

Fig. ^ shows the energy, £(0), and the density D{0, 0) for bottom states as a function of h. It is seen from the figure 
that for h ^ 2 the density falls well below its critical value 0.05, which supports our above conjecture, Sec. |l^, that 
the AEM provides reliable results for h ^ 2. With increasing h the angular momentum in the bottom state, L, also 
increases. It is a striking feature of the data that only states with L = 3n, n > 2, reach the spectrum bottom (we 
cannot make definite conclusions about the L = 3 state since it does nor reach the bottom for a = — 1/3 anyons, 
Fig. |l|). The bottom states described by L-even and L-odd polynomials alternate. We attribute the periodicity in 
L to the superselection rule originating from the combination of the space and permutation symmetry. Indeed, we 
observe this periodicity in the semiclassical region where the quantization rule includes-angular integration between 
two exchange points separated by the angle 27r/3 rather than the usual 27r integration.E^I 

Polynomials Pl.m with L ~ 6M play a special role in the class of L-even polynomials. All of them obey the hard 

L /3 

core constraint. Indeed, the polynomial Pgm.m (^12, 2^23, 231) vanishes as Zj/ each time when one of its arguments, 
Zji, turns into zero. From the standpoint of the general theory,0 hard-core functions are the only "legitimate" wave 
functions of an anyon system. The exponent L/3 is the maximum order of the zero for a wave function with the 
angular momentum L, and this maximum is achieved only for PeM,M polynomials. Therefore, Coulomb repulsion is 
strongly suppressed for these polynomials. Polynomials WPqm,m play a similar role in the class of L-odd polynomials. 
It is a remarkable fact, that in the h ^ 2 region all bottom states are either PeM,M or WPqm,m polynomials. 

Using the polynomials Pq^m, M = and 1, as an example, we can follow the way in which hard core states move 
to the spectrum bottom when h increases. For small h values the low energy component of the L = 6 doublet has a 
pronounced maximum near r = as it is seen in Figs. ^ and ||. For h < 2 the low-energy component is close to ^'5,0; 
for h = the overlap is 0.96. For h 2 wave functions of both components are strongly mixed, and they show similar 
distributions of the density. Fig. ^. For h > 2 the function ^Pg,! wins the competition. For h — 3 it dominates in the 
low-energy state; the overlap is 0.97. The density DQM,M{'r,0) has a single maximum for each value of M. One can 
obtain a simple analytic expression for the position of the maximum by averaging the density over r^. This latter 
function, Dl.m{'>', 0), reaches the maximum at — V2L, and the maximum of D^^Mif, 0) is very close to this value. 
More detailed information on the nature of the k = bottom states comes from the correlation functions which are 
discussed in Sec. ^ below. 

For comparison, in Fig. |^ is also shown the h dependence of the energy of a conventional magnetoexciton, emo(^), 
with the momentum k = 0. In the limit /i — > this energy exactly coincides with the energy of an exciton in the 
many-electron system because of the hidden symmetry inherent in the problem; see Refs. ^J[and |l^ and references 
therein. In the region of ft, ^ 1 the accuracy of the AEM is low. However, it increases for h ^ 2 when the charge 
fractionalization becomes important. In this region £mc(^) follows the usual Coulomb law, £mc(^) ~ whereas 
for many-electron systems the dependence of the exciton binding energy on h also is close to a Coulomb law, but the 
numerator is considerably less than 1 since the electron density distribution has a width about h. The magnetoexciton 
and AEM approaches are exact in the opposite limits. The results should be matched in the intermediate region at 
/i « 2. 

Above in this section we discussed the electron density distribution only for k = anyon excitons. The charge 
fractionalization manifests itself for these excitons in the crater shape of D\{r,Q). However, the most spectacular 
manifestation of the charge fractionalization can be expected in the large k region, fc ^ 1. Indeed, the exciton dipole 
moment d differs from k only by the rotation by 7r/2, Sec. |l[ Therefore, one can expect that with increasing k the 
electron density splits into bundles, their charges being multiples of 1/3. The splittpg of the electron shell into two 
well separated quasiparticles has been observed previously for a two-semion exciton.ll2l For a three-anyon exciton the 
pattern are much more impressive. For the bottom state, they are shown in Fig. ^ for ft, = 3 when the criterion of the 
large electron-hole separation is fulfilled. The distribution which is cylindrically symmetric for k — Q transforms with 
increasing k into a distribution with a single split-off anyon [k = 2 and 3). Two anyons constituting the exciton core 
show a slight but distinct splitting in a perpendicular direction. This core can be considered as an anyon ion. The 
core changes its shape with k but remains stable in a wide range of k. Finally, for rather large k values, it splits in 
the d direction as it is seen in the last figure, k ~ Q. The asymmetric density distribution for k 7^ arises completely 
due to the admixture of L-odd polynomials to the L = 6 state. 

The woil|Outlined profiles of the electron density seen in Fig. |^ may be smeared by the oscillatory screening inherent 
in IQL's.c2l Nevertheless, the basic pattern of the charge separation in an exciton should strongly influence the k 
dependence of the magnetoroton-assisted recombination processes since charge density excitations are left in a crystal 
afterwards. 
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VI. PAIR CORRELATION FUNCTIONS 



One can see in Fig. || that L — 3 state is the first state which shows a crater-hke shape of the density D(r, 0). This 
shape indicates the existence of the anyon substructure of an exciton as it was argued in Sec. 0. It is typical of all 
bottom states with L > 6. In this section we compare properties of V^s^o and '^q i states and show that despite the 
similarity in the shape of the density, they differ critically in the shape of the radial pair correlation function w{r). 

It is convenient to use the square of the wave function of Eq. (||) averaged over the hole coordinate. Using Eqs. (|^) 
and and performing Gaussian integration over Vh, one can write the following equation for the averaged '^qm.m 
function: 

|*6Af,Af (^12, ^23, 23i)P oc (riarzargi)-*^'^ expi-^irj^ + + rji)} ■ (41) 

Choosing ri2, ris, and ■0, where ip is the angle between the vectors ri2 and ri3, as independent variables, rewriting 
|^6M,Mp in terms of these variables, and looking for the maximum first over cosip and then over ri2 and ri3, one 
finds that the maximum of |5'6M,mP is reached for a configuration of an equilateral triangle with 



„2 



' 23 



= rl, ^ 6L . (42) 



In the semiclassical limit, L 3> 1, this result coincides with Eq. (|3q). 

To find the most probable configuration for the ^'3^0 state it is convenient to work in ^ variables. Sec. |l| and 
perform averaging over ■ Simple calculation shows that 



1*3,0(^12, ^23, ^3i)P (X ei&e^ exp{--(e? + + el)} ■ (43) 



The maximum of this expression under the constraint of Eq. (^ can be found in the same way as for Eq. (|l|). Finally, 
one recovers an equilateral triangle configuration with r^j — ''23 — ''31 = 18. This result coincides with Eq. (^2|) for 
i = 3. 

Because arbitrary L-odd polynomials have the form Pl,m — WPl-3,m, Eq- (p^), and the equilateral triangle 
configuration is optimal for each of the multipliers, it is optimal also for their product, Pl.m- Therefore, the most 
probable configuration has the same shape of an equilateral triangle both for L-even and i-odd states. 

To reveal a striking difference in the properties of '^q^i and ^'3^0 states, one can calculate the radial pair correlation 
functions w(r): 

w{ri2) = j |^'(R.,P, {zjiW ^(ri2 + r23 + rgi) dR dp dv2z dr^n . (44) 
Substituting j^'p from Eq. (||) results in 

WLM{ri2) ^ J^^J y c?i"23 dr3i|Pi,Mpexp{-^^r|, +iq- (ri2 +r23 + r3i)} . (45) 

Straightforward calculation shows that 

^6,i(0 = {1 + exp(-rVl2) , (46) 



^3.0(0 = {1 + 6^'/^^'^ exp(-rVl2) . (47) 

Both functions are shown in Fig. |^. The function we,i{r) has a hard core behavior, whereas W3fi{r) does not vanish 
at r = 0. On the contrary, it reaches its absolute maximum at this point. The second maximum is by the factor 
0.97 lower than the main one. Therefore, function q violates the hard core constraint. Of course, all bottom-state 
L-odd polynomials with L > 9 show a hard core behavior. 

The above conclusion on the triangular shape of the most probable configuration of the ^'3^0 exciton implies existence 
of the 120° correlations in the density correlation function. 

To investigate these correlations it is convenient to work in the variables. The two-particle correlation function 
W3,o{(p) depending on the angle (p between vectors and ^2 can be written as: 
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cx / dad6(Ci6)'(e? + e2 + 2a6cos(p)2 

Jo 

X exp{-(C? + el + 66 cos</p)/3} . (48) 

The delta function takes into account the constraint of Eq. (^, and ^^^2 standing for the angle between and 
^2 expressed in terms of the coordinates of these vectors. Integration over the variables in Eq. (48) results in the 



averaged correlation function, the main contribution comes from the area ~ Cl ~ ~ 6 where |^3,oP reaches the 
maximum. The last integral can be performed in polar coordinates, £,i — ^ cos(6'/2) , 6 = ^ sin(0/2) , < 6 < tt , 
and the normalized function w^^oif) takes the form: 

/ N 81 r sin^ 61(1 + sine cos (^) 

W3oi^) ^7—1 de — , . ^ . (49) 

Stt Jo {2 + sin 6 cos ip)^ ^ ' 

The integration over 6 can be performed analytically, however, the final expression is rather cumbersome. The result 
is plotted in Fig. ^. Two distinct maxima in the vicinity of 2tt/3 and 47r/3 reveal 120° correlations in the V^s^o state. 
The maxima are shifted from the angle 27r/3 since the function ^'3 vanishes for the coUinear configuration of anyons, 

ll = -|2, ^3 = 0. |-. 

Apalkov and Rashballil have found the density-density correlation function, uifg (</?), for the = 2/3 IQL with a single 
extra electron. Calculations were performed in the spherical geometry for the (iQp)max ^ Lqp ~ 3, {Lqp)z — Lqp 
quantum state. The quantity (iQp)max — Lqp, which is the difference between the three-quasiparticle angular 
momentum in the spherical geome.try, Lqp, and the maximum value of this momentum, should be compared to the 
exciton angular momentum LJIj'EZI It is convenient to introduce mean values 

/3l(l?) = j \^{ljJ,ljJ2,---i^N)\'^ dljJ2---dljJN (50) 

and 

P2(t^,'^ ) = / |*(u;,u; ,0^3.. .ljn)\^ du}3 . . . dujN ■ (51) 



Here Wj(i?, (p) are unit vectors designating the positions of the electrons on the sphere. Electron density depends only 
on the polar angle and equals ni{'d) = Npil-d), where N is the number of electrons. If one introduces the deviation, 
An(u)) = n(u)) — ni('d), of the density from its mean value, the density-density correlation function can be written as 

wi%{{>,(p-ip')^ < An{d,(p) An{d,(p') > ^ ^N(N ~ 1) p2{'&, ip ~ p ) - nji-d) . (52) 



In Eq. ( |52| ) the polar angles of the vectors (jJ and u; are chosen equal, ^ -d . Therefore, the correlation function, 
ifg g(i9, (p), is the function of the azimuth (p and depends on as a parameter. In Eq. ( |5^ ) the singular term proportional 

to 5{lj — u; ) is omitted, as usual, since it makes no contribution to g for uj ^ uj . The results of calculations are 
presented in Fig. ^ for three values of iS. The value of ?? = 0.8 is close to the maximum of fii (i?). 

Since the anyon exciton wave function 'i'3.0 does not depend on h, the data of Figs. || and g can be compared. There 
is a striking similarity between them. All three curves in Fig. ^ show flat but distinct maxima near the same values 
of the argument, ip = 2w/3 and 47r/3. These maxima are much flatter than the maxima in Fig. This difference can 
be attributed to the smearing originating from the quasielectron form-factors. Nevertheless, there is no doubt that it 
is the three quasielectron substructure of the many electron state L — 3 which manifests itself in 120° correlations. 

Therefore, the state \I'3^o shows properties intermediate between an anyon exciton and tight exciton. Unfortunately, 
reliable calculation of 'w'^'^{'d, (p — ip ) can not be performed for the (iQ£;)max — Lqe — 6 states since their size exceeds 
the accessible sphere size. 



VII. COMPARISON WITH FINITE SIZE CALCULATIONS 



In this section we make a comparison of the AEM with computational results for excitons in the spherical geometry. 
Three-anyon excitons are expected to appear in two charge conjugate IQL's, ly = 1/3 and i/ = 2/3. Their properties 
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in the semiclassical region, h ^ I and L ^ 1, should be identical. For intermediate h and L values, excitons,pf 
these two liquids are expected to show rather different properties. The data obtained by finite-size computationalj'EZI 
substantiate these expectations. 

There are different reasons why — 1/3 and 1^ — 2/3 excitons are expected to have different properties in the 
intermediate region. A simple electrostatic reason is a different nonlinear screening of the unit positive charge by 
these two IQL's. The effect of the screening on the exciton energy spectrum was discussed in Ref. However, there 
are also different mechanisms resulting in the difference in the properties of these excitons. They are related to rather 
different form- factors and different statistical properties of the quasielectrons of the = 1/3 and i/ = 2/3 IQL's. 
Neglecting spin effects, we can consider electrons as spin-less fermions. Then, because of the charge symmetry, the 
properties of the quasielectrons of the v = 2/3 IQL are identical to the properties of the quasiholes of the 1^—1/3 
IQL. Therefore, we will start with a comparative study of quasielectrons and quasiholes of the v — 1/3 IQL. 

The form-factors of charged particles of the v — 1/3 IQL were investigated in a number of papers. The data are 
summarized in Refs. |3^ - 13^. It is known that a quasihole has a narrow profile with a radius up to two magnetic lengths. 
The density decreases away from the center of a quasihole nearly monotonically. On the contrary, a quasielectron has 
a pronounced density dip at the center, the density maximum at about two magnetic lengths, and the radius of about 
four magnetic lengths. Therefore, the model of point anyons developed above matches much better the excitons of 
the u ~ 2/3 IQL than the excitons of the ly — 1/3 IQL. 

Another aspect of the problem is related to statistical properties of anyons. In the AEM the effect of the fractional 
permutational statistics of anyons was taken into account by including the factor Ilj/l^i)" i^^o the Halpcrin pseudo- 
wave-function '^LM J Eq. (^). It was shown in Sec. |^ that the effect of this factor can be neglected for h > 2, i.e., 
in the region where the AEM is expected to be valid. However, one should also take ipto account the nontrivial 
populational statistics of anyons. This can be done using the theory of co j gip osite fermional3 and the approach to the 
dimensionality of the quasiparticle space based on exclusion statistics.EJC3 Bosonic Haldane dimension, dgp, is an 
effective number of the single-quasiparticle states defined in such a way that the usual Bose distribution 

WB{d^P,NQp) = + Nqp - l)!/(d|p - ly.Ngpl 

results in the correct number of states in the Hilbert space of Nqp quasiparticles, W{Nqp) — Wpidgp, Nqp). 
The number of quasipaxticle states, W{Nqp), can be found by counting the number of states in the Hilbert space 
of composite fermions. If quasiparticles obey Bose statistics, dgp does not depend on Nqp. By counting the 
composite fermion number-of-states in the spherical geometry, it was shown in Ref. ^ that for ly = 1/m, m is an 
integer, the bosonic dimension of quasiholes equals 



d^H = N+l, (53) 



and the bosonic dimension of quasielectrons equals 

d^E = {N+l)-2{NQE-l) . (54) 

Here N is the number of electrons. These equations are consistent with the diagonal coefficients of [the_exclusion 
statistics of quasiholes and quasielectrons {1/m and 2 — 1/m, respectively), found by different authors.c3c3 

Equations ( |5^ ) and (jsj) indicate that quasiholes show a bosonic behavior, while quasielectrons are subject of 
the constraint which is even more restrictive than the constraint imposed by Fermi statistics. These conclusions 
were supportedE-3 by counting the maximum value of the angular momentum, (iQp)max , for the system of Nqp 
quasiparticles. It is equal to 

(iQH)max = NqhN/2 (55) 

for Nqh quasiholes, and to 

{LqeU.. = NqeN/2 - NqeINqe - 1) (56) 

for Nqe quasielectrons. Since N/2 is the angular momentum of a single quasiparticle j^^ Eq. ( |55| ) confirms the 
bosonic behavior of quasiholes, whereas Eq. -confirms the existence of the restriction on the population of single- 
quasiparticle states by quasielectrons. He et alSIl were the first to discover this restriction by means of numerical 
calculations and to propose Eqs. (^3|) and (|5^). They attributed the restriction to the hard core constraint for 
quasielectrons having a dynamical (short-range repulsion) rather than statistical origin. 



For a macroscopic system, iV 3> 1, equations (53) and ( |54| ) give coinciding results, dg^ w w N, in the dilute 
limit, Nqp <^ N. However, in an exciton the quasiparticles are confined inside the volume about 7rr^, where r is 
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the exciton radius. Therefore, the second term of Eq. ( p4[ ) which signifies the deviation from the bosonic behavior of 
quasielectrons, can be of importance. 

Quasielectrons of the = 2/3 IQL are described by Eqs. ( |53| ) and (|55|). Therefore, the population of single- 
quasiparticle states obeys the Bose statistics, and Eq. (y|) for wave functions is absolutely adequate since it includes 
polynomials Pl,m symmetrical in variables Zj. For this reason, and also taking into account narrow form-factors of 
quasielectrons of the v — 2/3 IQL, one can conclude that the ly = 2/3 IQL is the best candidate for comparison with 
the AEM. We believe, that the criterion h^2, established in Sec. ^ by evaluating the density D\{r, 0), is applicable 
to the 1/ = 2/3 IQL. 

The situation is more involved for the J/ = 1/3 IQL. Quasielectrons in this case are described by Eqs. ( p^ and ([5^). 
One can apply Eq. ( p4| ) to the area inside an exciton and order that dq^ > Nqe = 3, which results in > 6. The 

number of electrons inside an exciton can be evaluated as iV ~ 7rr^/27r. Evaluating the exciton radius r as r « ri2/V2> 
and using the classical equation (p5|), one comes to the criterion h> 5. A wide quasielectron form- factor imposes a 
similar restriction on h. Since the reliable finite-size computations can be performed only for /i 2, the prospects for 
a quantitative comparison of the results obtained by both approaches seem less favorable for the v ~ 1/5 IQL than for 
the V — 2/3 IQL. However, we feel that the above rigid criterion relaxes considerably when a qualitative description 
of the spectrum-bottom states is concerned. 

We are now in position to compare the basic results of the AEM with the computational results for finite-size 
systems. The basic statement of th e . th eory of anyon excitons, that the charge fractionalization results in a multiple- 
branch energy spectra of excitons ,tilE3 was confirmed by finite-size computations performed for both = 1/3 and 
V — 2/3 IQL's. Actually, it provides the basic idea for representing the energy spectra obtained in the spherical 
geometry for a discrete set of L values in the form of exciton branches £\ (k) . 

We begin with summarizing some results obtained in Ref. |l^ by finite-size calculations for excitons of the = 2/3 
IQL and compare them with the results obtained in this paper in the framework of the AEM. Classification of excitons 
in terms of tight and anyon excitons is used, cf. Sec. |. 

1. Number (^.expiton species. Anyon exciton is a bound state of a three-quasielectron complex from the low-energy 
(anyon) sectorEZlEj and a hole. The angular momentum of the exciton, L, is equal to L* = {LQE)m&x — Lqe, where 
Lqe is the angular momentum of the complex, and (iQ£;)max is the maximum value of this momentum which can 
be found from Eq. (55). Therefore, the number of exciton species is equal to the number of the three-quasielectron 
complexes in the low-energy sector, and L* should be compared to the exciton angular momentum L of the AEM. 
Only a single three-quasielectron complex exists if the angular momentum L* equals L* — 0, 2, 3, 4, 5, or 7. The 
complex L* = 1 is absent because composite fermions obey Fermi statistics. For L* > 6 the number of L-even states 
increases by one (see Fig. 1 in Ref. p^ . 

These properties are in a complete agreement with the classification of the polynomials Pl,m of Sec. ||, Eqs. (H) 
and ( p^ ) . The L = 1 exciton is absent because of the constraint of Eq. (H) . 

2. L = branch. Despite the fact that the L = exciton appears, according to its quantum number, as the first 
state in the series of anyon excitons, it possesses rather special properties. This exciton originates from the L* = 
three-quasielectron complex which is quite dense, and the energy of this complex is high. It is nearly the same as the 
energy of some states from the next sector. These data imply that the L — anyon exciton actually merges with the 
L — tight exciton. Because of these arguments, the L = exciton was assigned in Ref. |l^ as a tight exciton rather 
than an anyon one. This assignment is supported by an independent argument. L — excitons of the v = 2/3 and 
1/ = 1/3 IQL's show nearly identical properties, whereas the anyon exciton assignment of the latter entity is excluded 
by the symmetry arguments based on the composite fermion theory.E3 

The shape of the L = exciton density distribution of Fig. ^ is in agreement with this assignment. The magnitude 
of Do{0, 0) exceeds the maximum density compatible with the AEM; see the discussion at the beginning of Sec. 0. 

3. Electron density. For a system with a single extra electron against a background of the ^' = 2/3 IQL, the electron 
density has a pronounced maximum at r = if L* = 0, 2, or 4. The L* = 3 state is the first state with a crater-like 
electron density distribution having a density dip at r = (see Fig. 2 in Ref. |l^). The crater shape of the density 
indicates the charge fractionalization. 

These data are in a qualitative agreement with the electron density distributions of Fig. ^ for excitons having 
angular momenta L — 0, 2, 3, and 4. The general shape of the curves is the same, but there are differences in the 
magnitudes of the densities at r = 0. Since wave functions of these excitons do not depend on /i, electron densities 
for the exciton and three-quasielectron states with the angular momenta L = L*, respectively, can be compared. 

4. Bottom states. L = exciton remains the bottom state in the whole area h ^ 2. L — 2, L — 3, and L — 6 
excitons move down to the spectrum bottom, but they start competing with the L — exciton only for h « 2.5 when 
the accuracy of finite-size calculations becomes aipbiguous. L = A exciton is also seen in the low-energy part of the 
spectrum but never reaches the spectrum bottomO^EJ 

These data are in agreement with Fig. where the sequence of the first bottom states includes L — 0, 2, ,3 and 6. 
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5. Charge fractionalization: density correlation function. It was shown in Sec. Vl that W3fi{ip) reveals 120° corre- 
lations both for the many-electron and AEM wave functions. These correlations signal the charge fractionalization. 

6. Intrinsic angular momenta of anyons. There is a convincing one-to-one correspondence between the excitons 
of many electron systems having angular momenta L < 6 in the spherical geometry and the excitons of the AEM 
with the same values of the z projections of the angulax-momcntum. Therefore, our data provide no indication of the 
existence of the intrinsic angular momenta of anyons.EJ On-the contrary, our data are in agreement with the recent 
conjecture on the absence of anyon spins in the plane limitEd |— in 

The AEM predicts identical exciton spectra for the v = 2/3 and v — l/2> IQL's. However, finite-size calculational3'EZl 
result in a rather different symmetry of the low-energy exciton states for ft, < 2. The exact classification of the exciton 
states based |OIl the composite fermion theory shows that only L = 3 and L > 5 anyon excitons can exist in the 
I' = 1/3 IQL£3 This conclusion is supported by numerical data. Therefore, a challenging question arises: Why L = 2 
and L — A excitons of the AEM are missing from the many-electron spectra of the = 1/3 IQL? We argue, that these 
excitons are excluded because of their small bosonic dimension which cannot accommodate three quasielectrons. 

In what follows we compare the results of the finite-size computations of Ref. [l^ for the = 1/3 IQL with the AEM 
data. 

1. Excluded states. Finite-size calculations and composite fermion theory show that anyon excitons of the = 1/3 
IQL can only posses angular momenta L = Z and L > 5. All excitons with L = 0, 1, 2, and 4 can only appear as 
tight excitons. L — {) tight exciton forms the spectrum bottom for h ^ 2, but there is no low-energy L — 2 and L = 4 
excitons. 

It seems probable that there can exist only one, L — 0, tight exciton near the spectrum bottom. Therefore, it is 
necessary to understand why L = 2 and L = 4 excitons of the AEM, which possess low energies, do not appear as 
anyon excitons in many-electron systems. The criterion dq^ > Nqe results in the minimum anyon-exciton radius of 

^min ~ 2^3. It is seen in Fig. |^ that the criterion r ^ rmin is violated for L — 2 and L — A excitons. Therefore, we 
attribute the exclusion of the L — 2 and L = 4 exciton states from the many-electron spectrum to the reduction of 
the bosonic dimension of the quasielectron space because of the second term of Eq. (p^). 

2. Bottom states. Tight L = exciton remains the bottom state up to ft « 2 when the L — 3 anyon exciton reaches 
the bottom. L = 5 and L = 6 excitons start to compete with the L — exciton only for ft w 3 when the accuracy of 
finite-size calculations is low. 

The data are in agreement with Fig. ^. 

3. Electron density distribution for k = excitons. Tight excitons have narrow electron density distributions. The 
density of the L = exciton is sharply peaked at r = 0. L — i and L — 5 anyon excitons show pronounced density 
dips at r = and two maxima of the density. It was hypothesizedtj that the r = dip originates from the dip in the 
quasielectron form-factor, whereas the two-maxima shape indicates the existence of a two-anyon core and one split-off 
anyon. It is remarkable that the density distribution in the L — 3 and L = 5 excitons only weakly depends on ft. 

Weak ft dependence of the density distribution for L ~ 3 and L = 5 excitons is in agreement with the existence of 
a single polynomial Pl,m for L < 6, Sec. ||. Splitting of the exciton shell of an k = exciton can not be described 
within the framework of the point-anyon AEM. 

4. Dependence of the density distribution on k. When k increases, the r — dip in the electron density of the 
L = 3 cxcitotujtransforms into a narrow maximum, and the density distribution acquires a three maxima shape. It 
was proposccEj that in this region of k values an exciton consists of a single-anyon core and a two-anyon shell. 

Splitting off of anyons from the core with increasing fc is in a qualitative agreement with Fig. ^ 



VIII. CONCLUSION 



The model of anyon exciton developed and solved in the paper includes three quasielectrons (anyons) and one hole. 
It is applicable to exciton spectra of two charge conjugate IQL's, v = 1/3 and v — 2/3, and is exact in the limit of a 
large separation between electron and hole confinement planes, h ^ I. 

Anyon excitons possess multiple branch energy spectra, ex{'k). An exciton is described by a 2D momentum k and 
two internal quantum numbers A which numerate exciton branches. One of these quantum numbers can be chosen as 
the exciton angular momentum, L, in the k = state. 

A full set of basis functions was chosen with a proper account of the magnetic translational symmetry and permu- 
tational symmetry. The functions include two types of polynomials symmetric in anyon permutations, one of them, 
apparently, was considered for the first time. Analytic expressions for all matrix elements were derived. As a result, 
exact expressions for the energy spectrum of a four-particle system were found for k = 0. All exciton states with 
even angular momenta L < 6 and odd momenta L < 9 are non-degenerate, and their wave functions are completely 
determined by symmetry requirements. This property manifests itself in finite-size data in a weak ft dependence of the 



15 



electron density distribution. Properties of k 7^ excitons were investigated by numerical solution of the Schroedinger 
equation. 

The AEM is not applicable for small h values, h < 21. However, the analysis of the internal criteria of the AEM 
and comparison with the finite-size data show that it gives satisfactory results for /i « 2Z, especially for excitons of 
the v = 2/3 IQL. Finite-size computations which are reliable for h < 21, and the AEM, whose accuracy increases 
with h, constitute two supplementary approaches which match at « 21. According to the AEM, exciton states with 
the angular momenta obeying the superselection rule L — 3m, where m > 2 is an integer, form the sequence of the 
bottom states for h > 21. All these states are hard core states, i.e., their wave functions turn into zero if any two of 
the anyon coordinates coincide. The equilateral anyon configurations are the most probable ones. Angular momenta 
L of the bottom states increase as h'^, and the size of these states as h. The tight L = state is the bottom state for 
h < 21, whereas L = 3 state, which possesses the properties of both anyon and tight excitons, can reach the bottom 
for h w 21. 

Anyon substructure of excitons manifests itself in their different properties. First, the multiple branch energy 
spectrum originates because of the existence of the internal degrees of freedom. Second, for k = excitons a crater- 
like shape of the electron density distribution with a pronounced dip near the hole signals the existence of the charge 
fractionalization. Third, 120° correlations in the electron density were found both in the framework of the AEM and 
in finite-size data. Fourth, for k 7^ excitons a spectacular splitting of the electron density into bunches provides a 
direct manifestation of the fractional charge substructure. 

AEM results in identical low-energy spectra for the excitons of the h' — 1/3 and h' — 2/3 IQL's. This result 
is definitely correct in the h 3> I limit. However, finite-size data suggest that in the intermediate region, h ~ I, 
the v = 1/3 IQL possesses less anyon exciton branches than the ly = 2/3 IQL. Following Ref. 0, we attribute this 
phenomenon to the difference in the populational statistics of the quasielectrons of the charge conjugate IQL's. Anyon 
statistics enters into the theory in two ways. First, through the fractional-power factors in the wave function; they 
ensure correct interchanging statistics. Second, through the populational statistics. We have shown that the first 
mechanism can be neglected for h ^ 21, whereas the second one is of a critical importance in the intermediate region. It 
follows from the composite fermion theory that quasielectrons of the 1^ — 2/3 IQL obey bosonic populational statistics. 
As a result, quantum numbers of anyon excitons found from the AEM and from the finite-size data exactly coincide. 
On the contrary, the bosonic dimension, dg^, of the quasielectron space of the v — 1/3 IQL rapidly decreases with the 
number of quasielectrons, Nqe- The condition dg^ > Nqe, written for the area about the exciton size, eliminates 
the exciton states of the AEM with the angular momenta L = 2 and L = 4 from the many-electron spectrum and 
brings in agreement the AEM and many-electron data. Therefore, for h k, 21 the AEM supplemented with anyon 
statistics arguments matches the many-electron data for both IQL's and sheds light on the origin of the difference in 
their exciton spectra. The above arguments are rather general and can be applied to different problems of the theory 
of the bound states of several anyons. 

Comparison of the AEM and finite-size data for v — 2/3 IQL provides no indication of the existence of the intrinsic 
angular momenta (spins) of anyons. 

In conclusion, the anyon exciton model unveils the general pattern of the exciton spectra of IQL's. It predicts the 
properties of excitons in the large h limit, I, and is in agreement with finite-size data in the intermediate region, 
h fa 21. The theory suggests that the region /i > 2/ is most favorable for investigating the anyon substructure of 
excitons in optical experiments. 
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FIG. 1. Ground state energy, e(k = 0), plotted versus separation, h, between electron and hole confinement planes. A and 
B - statistical charges a — —1/3 and a = 0, respectively. A - triangles show the points where the angular momentum L changes 
from zero to 2, and then to 6 and 9. For the quantum states separated by full dots on the curve B see Fig. C - classical limit. 
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FIG. 2. Anyon cxciton dispersion law e{k) for two values of h. For h = 2, the negative dispersion arises because of the 
mutual rcpulsiou of L = 2 and L = 3 branches. Level splitting near avoided crossings becomes tiny with increasing h. Numbers 
show L values, h in units of I. 
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FIG. 3. Axisymmctric electron density distributions Z)L(r, 0) for an k = cxciton for the states with L < 6. Two L 
states are shown for h = 0; the density distribution in the lower state is shown by a solid line. Numbers show L values. 
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FIG. 4. The energy £(0) and the electron density D{Q, 0) at the point Vh where the hole resides plotted vs h for the ground 
state of an exciton with A: = 0. The ground state energy of an anyon exciton is shown by a solid line; the dots on it show the 
positions of the intersections between the energy levels with different L values. For comparison the energy of a conventional 

magnctoexciton e-mcih) with fc = is shown by a dashed line. Numbers near the D{0, 0) curve show the L values. Only the 
states with L = 3m reach the spectrum bottom (as an exclusion the state L = 2 appears as a bottom state in an extremely 
narrow region of the h values) . 
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FIG. 5. Electron density distribution D{r, 0) for the L = 6 states with A; = for three values of h. The density in the lower 
energy state is shown by a solid line. Consecutive numbers, A, of the energy levels are shown near the curves. 
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FIG. 6. Electron density distribution in an ariyori cxciton for different values of k. A hole is at the origin, the x axis is 
chosen in the d direction. The center of the electron density distribution is at a; = k, y=0. The data were obtained with the 
44 polynomial basis, L < 21; data for the 70 polynomial basis, L < 27, also show the change in the shape of the exciton core 
for large k values, x, y and h in units of I; k in units of density in units of 
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FIG. 7. Radial anyon pair correlation function w{r) for the states ^'e.i and ^s^; k = 0. 
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FIG. 8. Normalized anyon pair correlation function W3fi{(p) found in the anyon exciton model; k = 0. 
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FIG. 9. Electron density correlation function wg'^Q{-d, ip) for the v — 2/3 IQL with a single extra electron (A'^ = 15, the flux 
2S — 21). The correlation function is plotted as a function of azimuth (p for three values of the polar angle = 0.6, 0.8, and 
1.0. 1? — 0.8 is close to the maximum of the electron depsity. Correlations between quasielectrons manifest themselves in the 
maxima near (f = 2tt/3 and Att/3. (Apalkov and RashbaHl) 
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